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An efficient algorithm for classification of binary self-dual codes is presented. As an 
. application, a complete classification of the self-dual codes of length 38 is given. 
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^ ! 1 Introduction 

The self-dual codes form one of the important classes of linear codes because of their rich 
algebraic structure and their close connections with other combinatorial configurations like 
block designs, lattices, graphs, etc. 

The classification of self-dual codes began in the seventies in the work of Vera Pless 
[20] . where she classified the binary self-dual codes of length n < 20. The method used in the 
beginning remained essentially the same throughout the succeeding classifications. This is a 
recursive classification which proceeds from smaller to larger length and codes are classified 
up to equivalence. The process begins with the formula for the number of all self-dual codes 
of length n called a mass formula. The number of the self-dual binary codes of even length 

nis N(n) = nrii _1 ( 2i + 

Throughout this paper all codes are assumed to be binary. Two binary codes are 
called equivalent if one can be obtained from the other by a permutation of coordinates. The 
permutation a € S n is an automorphism of C, if C = o~(C) and the set of all automorphisms 
of C forms a group called the automorphism group of C, which is denoted by Aut(C) in this 
paper. If C has length n, then the number of codes equivalent to C is n!/|Aut(C)|. To classify 
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self-dual codes of length n, it is necessary to find inequivalent self-dual codes C\,...,C r so 
that the following mass formula holds: 



In the survey [T2J Huffman summarized the classification of all binary self-dual codes 
of length n < 36. As the complete classifications for lengths 34 and 36 was not done at that 
time, we present here a new version of Table 1 from |12j . but instead of the question marks 
we put the correct numbers. Moreover, we extend the table with two more lengths - 38 and 
40. The number of all inequivalent singly-even (Type I) and doubly-even (Type II) codes 
is denoted by (jj and jjjj, respectively. In the table d max j (d max ,ii) is the biggest possible 
minimum distance for which Type I (Type II) codes with the given length exist, and $max,l 
(Umox,//) is their number. 

The classification of the self-dual codes of length 34 was given in [3j. Using the self- 
dual [34, 17, 6] codes, Melchor and Gaborit classified the optimal self-dual codes of length 
36, namely the [36,18,8] codes [19]. Recently, Harada and Munemasa in |11| completed the 
classification of the self-dual codes of length 36 and created a database of self-dual codes [ID] . 
The doubly-even self-dual codes of length 40 were also classified by Betsumiya, Harada and 
Munemasa [2J. Moreover, using these codes, they classified the optimal self-dual [38,19,8] 
codes. A classification of extremal self-dual codes of length 38 was very recently obtained 
also in [1] by somewhat different techniques. 

Actually, the construction itself is easy, the equivalence test is the difficult part of the 
classification. Because of that, for the larger length the recursive constructions are used 
preferably heuristic, for building examples for codes with some properties. There are also 
many partial classifications, namely classifications of self-dual codes with special properties, 
for example self-dual codes invariant under a given permutation, or self-dual codes connected 
with combinatorial designs with given parameters. These classifications are not recursive but 
they use codes with smaller lengths that is why the full classification is very important in 
these cases, too. 

In this paper, we present an algorithm for generating binary self-dual codes which gives 
as output exactly one representative of every equivalence class. To develop this algorithm, 
we use an approach introduced by Brendan McKay known as isomorph-free exhaustive gen- 
eration [18] , The constructive part of the algorithm is not different from the other recursive 
constructions for self-dual codes, but to take only one representative of any equivalence class, 
we use a completely different manner. This approach changes extremely the speed of gener- 
ation of the inequivalent codes of a given length. Its special feature is that practically there 
is not equivalence test for the objects. 

As a result, the classification of all binary self-dual codes of length 38 is presented. The 
number of these codes is given in the following theorem and also listed in Table [TJ 

Theorem 1 There are 38 682 183 inequivalent self-dual codes of length 38. 




(1) 
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In Section [2] we present the theoretical foundations of our construction method. We de- 
scribe a recursive construction of binary self-dual codes. The aim is to obtain all inequivalent 
self-dual codes of dimension k on the base of the inequivalent self-dual codes of dimension 
k-1. 

In Section we describe the used algorithm. Codes which are equivalent belong to the 
same equivalence class. Every code can serve as a representative for its equivalence class. We 
use the concept for a canonical representative, selected on the base of some specific conditions. 
This canonical representative is intended to make easily a distinction between the equivalence 
classes. 

In Section S] we give the results with some more information about the obtained codes. 



Table 1: Binary self-dual codes of length n < 40 [12J 
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2 Recursive Constructions of Self-Dual Codes 

In this section we describe a recursive construction of binary self-dual codes. The aim is to 
obtain all inequivalent self-dual codes of length n = 2k and dimension k on the base of the 
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inequivalent self-dual codes of dimension k — 1. Theorem [3] gives us the possibility to develop 
such a recursive algorithm. Actually, this theorem presents a well known property of the 
self-dual codes but usually the authors state it in the case d > 4 (see for example [H]). Here 
we give a variant for all d but for that we need the following lemma. 

Lemma 2 If C is a binary self-dual code then C does not have more than two equal coor- 
dinates. In other words, if x^ = Xi 2 = ■ ■ ■ = Xi s for any codeword x = {x\,x%, . . . ,x n ) € C 
then s < 2. 

Proof. If n = 2 then C = i 2 = {00, 11} and C has only two coordinates. Let n > 4. Since 
the dimension of C is n/2 > 2, not all coordinates are equal and s < n. Suppose that s > 3, 
so Xi x = Xi 2 = Xi 3 for all codewords x € C. In such a case the vector y = {y\, ... ,y n ) € F?j 
of weight 2 with y^ = yi 2 = 1 will belong to C 1 - = C. In this way we obtain a codeword for 
which yi x = yi 2 = 1 but yi 3 = - a contradiction. Hence s < 2. □ 

Any self-dual [n, n/2, 2] code for n > 2 is decomposable as i 2 © C n _2 where C n _2 is a 
self-dual code of length n — 2. Furthermore, the number of the self-dual [n, n/2, 2] codes is 
equal to the number of all self-dual [n — 2, n/2 — 1] codes (n > 2). But even if a self-dual code 
of length n > 2 has minimum distance d = 2, it has two coordinates which are not equal. 

Theorem 3 Let C be a binary self-dual [n = 2k > 2,k,d] code and Co = {(x±, . . . ,x n ) £ 
C, x n _i = x n }. if the last two coordinates of C are not equal and C\ is the punctured code 
of Cq on the coordinate set {n — l,n} then C\ is a self-dual [n — 2, k — 1, d\ > d — 2] code. 

Proof. Let n > 2 and the last two coordinates of the code C are not equal. This means 
that there is a codeword y € C such that y n -i ^ Un- Therefore Co is a subcode of C with 
dimension k — 1, and (00... Oil) C. The punctured code of Co on the coordinate set 
{n — 1, n} is 



Since (00 ... 011) C, dimCi = dim Co = k — 1. Moreover, if x = (x%,X2, ■ ■ ■ , x ra _2, a, a), 

V = (yi, 2/2, • • • , 2/n-2, 6) G C then xij/i + x 2 2/2 H + x n - 2 y n ~2 + ab + ab = 0. Hence 

+ x 2 y 2 H h x n -2Vn~2 = for any two codewords (x%, . . . , x n _ 2 ), (yi, ■ ■ ■ , 2M-2) € Ci. 

It follows that the code C\ is self-orthogonal, and since its dimension is a half of its length, 



For our construction, we use the code C\ of dimension k — 1 to obtain a self-dual code 
of dimension k and length 2k. We describe this in the following corollary. Let C\ be a binary 
self-dual [2k — 2,k — 1] code with a generator matrix G\. Let us extend G\ with two equal 
columns (oj, 02, . . . , a^_i) T and consider the code Cq generated by the matrix 



Ci = {(xi,x 2 , . . • ,x n _ 2 ) I 2C = (x\,xv, ...,x n 



2,X n -l,X n ) G C }. 



it is self-dual. 



□ 
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We can choose the vector (a%, 02, ■ ■ ■ , dfe-i) such that (11... 11) G Co- Obviously, Co is a 
self-orthogonal [n = 2k, k — 1] code and its dual code is Cq = (Co, (00 . . . Oil), x) where 
x G Cq \ (C , (00 . . . Oil)) and wt(x) is even. Consider the last two coordinates of the vector 
x. If x n -\ = x n then (x±, X2, . ■ . , #0-2) £ Cj 1 = C\ and therefore x G (Co, (00 . . . Oil)). 
Hence x n _\ ^ x n . The next corollary follows immediately. 

Corollary 4 The code C = {C$,x) = Co U (x + Co) is a binary self-dual [2k, k] code. 

We will call the code C a child of C\, and Ci - a parent of C. Actually, we have two 
possible self-dual codes obtained in this way from C\, namely (Co, x) and (Co, (00 . . . 011)+x), 
but both codes are equivalent. 

Remark 1 To describe the search tree for our algorithm, we use the terms parent and child 
as it is usual in the literature on algorithms |15|, [T8] . These terms are used in a different way 
in some papers and chapters devoted to self-dual codes [13]. 

In this way all self-dual codes of length n = 2k can be constructed on the base only on 
the inequivalent self-dual codes of length n — 2. Indeed, if we take two equivalent self-dual 
codes Ci = C[ and use the same vector (01,02, ... , Ofc-i), we can obtain equivalent self-dual 
codes of length 2k via the described construction. Let G\ be a generator matrix of C\ and 
P be a permutation matrix such that G\P generates the code C[. Then the matrix 



generates a self-dual code equivalent to C, namely the code C. Furthermore, if ir G S n -2 is 
the permutation corresponding to the matrix P then tt G S n sends C to C where ir(i) = ir(i) 
for 1 < i < n — 2, and n (i) = i for i G {n — 1, n}. 

Let see now what happens if we take different vectors a and b of length k — 1 and use 
them in the described construction for the same self-dual [2k — 2, k — 1] code Ci with a gener- 
ator matrix G\. Consider the elements of the automorphism group Aut(Ci) as permutation 
matrices of order n — 2. To any permutation matrix P G Aut(Ci) we can correspond an 
invertible matrix Ap G GL(A: — 1,2) such that G' x = Gi-P = ApG\, since is another gen- 
erator matrix of C\. In this way we obtain a homomorphism / : Aut(Ci) — > GL(/c — 1, 2). 
Consider the action of Im(f) on the set Fg -1 defined by A(x) = Ax T for every x G F^ 1 . 

Theorem 5 [llj The matrices (Gi a T a T ) and (G\ b T b T ) generate equivalent codes if and 
only if the vectors a and b belong to the same orbit under the action of Im{f) on Fg -1 . 

Proof. Let the matrices (Gi a T a T ) and (Gi b T b T ) generate the codes Co and C , respectively, 
and a T = A P b T , where P G Aut(Ci). Then 
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Since Ap(G\ b T b T ) is another generator matrix of the code G , both codes are equivalent. 

Conversely, let Cq = G . It turns out that there is a matrix B G GL(k, 2) and an n x n 
permutation matrix P such that (G\ a T a T ) = B(G\ b T b T )P = BA P (Gi b T b T ). Hence 
BApGi = G\ and therefore BAp defines an automorphism of the code C\ = (G\). Since 
BA P b T , the vectors a and b belong to the same orbit under the action of Im(f) on 



F. 



fc-l 



□ 



Let now G 
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If P G Aut(Ci) and y = xP then 



be a generator matrix of the self-dual [n, n/2, d] code C. 
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o h 
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b T b T 



Hence the code C is equivalent to the code generated by the matrix G' 



y 1 

The following theorem reduces the number of the considered cases. It is a particular 
case of Theorem 1 from |16| . 



Theorem 6 If C is a binary self-dual [n = 2k > 2,k,d] code and the last two coordinates of 
C are not equal then C is equivalent to a code with a generator matrix in the form 



( x\... x k _ x 00 . 



G 



1 



\ 



h 



fc-l 



A 



\ 



(2) 



and the matrix {Ik-\\A) generates a self-dual [n — 2, k — 1] code. 



Proof. Take a generator matrix in systematic form for the code C\ obtained from C by the 
construction in Theorem [3j Then 



V 



xi 



xi 



'fc-l 



A 



Xk-l X k -1 



is a generator matrix of Cq in systematic form. Since Cq is a self-orthogonal [2k, k — 1] code, 
Cq = (Cq, (00 . . . Oil), y) where y G C \ Cq. Hence y n _i / y n and suppose y n -i = 1, y n = 0. 

Consider the vector x = (xi,X2, ■ ■ ■ , a?fc-i, 00 . . . 010). Obviously, x _L Go and hence 
x G Gq 1 . But x G , x (00... 011) + Cq, x (00... 011) + y + G and therefore 
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x G y + Co C C. It turns out that x € C and we can take a generator matrix for C in 
the needed form. □ 



If we start with a self-dual [2k — 2, k — 1] code Ci with a generator matrix in systematic 
form, for the construction in the last theorem we have to use a vector (xi,X2, ■ ■ ■ ,Xk—i) of 
odd weight. Then we are sure that the all-ones vector (11 . . . 11) belongs to the subcode Cq 
and the code C generated by the matrix ([2]) is a self-dual [2k, k] code. 

The last two theorems are very important for our search. We use them to write a 
recursive algorithm which gives us all self-dual [n,n/2] codes starting from the code %i. 

3 The Algorithm 

In this section we present an algorithm for generating self-dual codes of a given length n. We 
discuss the background and give preliminary definitions and notations. 

Let fifc be the set of all binary self-dual codes of dimension k (and length 2k). We 
consider the action of the symmetric group S2k on the set fifc, k = 1,2, ... . This action 
induces an equivalence relation in as two codes C\,C2 € Ofc are equivalent {C\ = C2) if 
they belong to the same orbit. Hence the equivalence classes for the defined relation are the 
orbits with respect to the action of the symmetric group. According to Theorem[3l if we take 
two nonproportional coordinates in any code belonging to k > 2, we can obtain a code 
from . Conversely, if we take all codes from Qk—i an d extend them using Corollary [5] in 
all possible ways, we will obtain all codes from Q^. The construction of the self-dual codes 
of dimension k using the codes from flk-i seems to be trivial, but actually this is a difficult 
problem and the question is how to find only the inequivalent codes. To do this, we develop 
a McKey type algorithm for isomorph-free generation |15} [T8] . 

The algorithm is an exhaustive search over the set of codes f2 = fii U 02 U ■ ■ ■ U fij 
and this set is our search space. The generation process is described by a rooted tree (or 
forest). The nodes in the tree are objects from the search space Q. Prom a self-dual code A 
of length 2k — 2 corresponding to the node A, we obtain self-dual codes of length 2k which 
are children of the code A. To construct the children, we use a generator matrix Ga of A 
in systematic form and the binary vectors from Fg 1 of odd weight. We denote the set of 
inequivalent children by Child(A). The elements of Child(A) correspond to the nodes of the 
next level which are connected to A by edges. To find only the inequivalent children, we use 
Theorem [5j Practically, the rule A — > Child(A) which juxtapose children to a code defines 
our search tree. The execution of the algorithm can be considered as walking the search 
tree and visiting all nodes through the edges. This can be done by a depth first search. We 
present a pseudocode of the algorithm in Table [2j 

The following definitions and theorems help us to explain the algorithm and to prove 
its correctness. 

Definition 1 A canonical representative map for the action of the group S2k on the set 
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is a function p : £2^ — > that satisfies the following two properties: 

1. for all X en k it holds that p(X) ^ X, 

2. for all X,Y <E Q k it holds that X ^ Y implies p(X) = p(Y). 

For a code C £ Clk, the code p(C) is the canonical form of C with respect to p. Analo- 
gously, C is in canonical form if p(C) = C. The code p{C) is the canonical representative of 
its equivalence class with respect to p. 

We can take for a canonical representative of one equivalence class a code which is 
more convenient for our purposes. Suppose that A is the canonical representative of this 
class with respect to a canonical representative map p. We can take B = A for a canonical 
representative for this class if we change the canonical representative map in the following 
way: p'(X) = B if X ^ B, p'{X) = p(X) if X ¥ B. According to Lemmata self-dual code 
does not have more than two equal coordinates, hence if k > 2 we can take for a canonical 
representative of any equivalence class a code for which the last two coordinates are not equal. 

Example 1 If we order the codewords in any code lexicographically and then compare the 
codes according to a lexicographical ordering of the columns, we will have one biggest code 
in any equivalent class. We can take this code as a canonical representative of its class. 

Definition 2 Let C be a self-dual code of length n and p(C) = C be its canonical form. A 
permutation (ftc & S n is called a canonical permutation of the code C, if <pc(C) = p(C). 

For a fixed canonical representative map p, the canonical permutation of C depends on 
the automorphism group of the code and the permutations from the coset (/>cAut(C) can also 
be canonical permutations since 4>ccr(C) = 4>c(o-{C)) = 0c(C) = p(C) for any a E Aut(C). 

We obtain the automorphism group and a canonical form of a given code A using a 
modification of the algorithm presented in the paper [5]. This algorithm gives the order of 
the group, a set of generating elements, and a canonical permutation. 

The parent test can be defined in the following way. Suppose that A is a self-dual code of 
dimension k — 1 and B € Child(A). This means that B is obtained from A by the construction 
of Theorem [6j The parent test for B depends on the automorphism group of the code B. Let 
ci and C2 be the coordinate positions of the code B for which 4>b{c\) = 2k,<f>B(c2) = 2k — 1 
where 0b is the canonical permutation of B. The corresponding two columns of B are not 
equal, because the last two coordinates of the canonical representative are not equal. We call 
the coordinates c\ and C2 special for the code B with respect to the canonical permutation 
4>b- If there is an automorphism a of B such that {cr(ci), c(c2)} = {2k — 1,2k} then the 
code B passes the parent test. In such case we can change the canonical permutation, taking 
4>bo~~ 1 instead of 4>b- Then the last two coordinates of B are special with respect to the new 
canonical permutation. So a code passes the parent test, if there is a canonical permutation 
<f> for this code such that the last two coordinates are special with respect to 4>. 

Obviously, the canonical representative B of one equivalent class passes the parent test. 
If B\ = B also passes the parent test, then there is a permutation <p : B — > B\ such that 
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1 < <p{i) < n — 2 for 1 < i < n — 2 and {<f>(n — 1), 4>(n)} = {n — 1, n}, where n is the length 
of the codes. This means that the parents of these two codes are equivalent, too. So we have 
the following lemma. 

Lemma 7 // B\ and B2 are two equivalent self-dual [2k, k] codes which pass the parent test, 
their parents are also equivalent. 



Table 2: The main algorithm 



Procedure Augmentation^: binary self-dual code); 
begin 

If the dimension of A is equal to k then 
begin 

U k := U k U {A}; 

PRINT (A,|JAut(A)); 
end; 

If the dimension of A is less than k then 
begin 

find the set Child(A) of all inequivalent children of A; 

( using already known Aut(A);) 
For all codes B from the set Child(A) do the following: 
if B passes the parent test then Augmentation^); 

end; 

end; 

Procedure Main; 

Input: U r - nonempty set of binary self-dual [2r, r] codes; 

k - dimension of the output codes (k > r); 
Output: Uk - set of [2k, k] binary self-dual codes; 
begin 

Uk '■= (the empty set); 

for all codes A from U r do the following: 

begin 

find the automorphism group of A; 
Augmentation^); 
end; 
end; 
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Lemma 8 Let A\ and A% be two equivalent self-dual codes of dimension r. Then for any 
child B\ of Ai which passes the parent test, there is a child B2 of Ai, equivalent to B\, such 
that B2 also passes the parent test. 

Proof. Let G\ be a generator matrix of A\ in systematic form and B\ be the code obtained 
from A\ and the vector a = (01,. . . ,a r ) by the construction described in Corollary H] and 
Theorem [6] Let B2 be the code generated by the matrix vr(Gi) and the same vector a, where 
7r G S2r and n(Ax) = A2. Obviously, vr(Gi) generates the code tt(Ai) = A2- Moreover, 
according to Theorem [5j B2 is equivalent to all codes obtained by A2 and the vectors from 
the orbit with representative a under the action of Aut(^2) on Fj. Since the codes B\ and B 2 
are equivalent, they have the same canonical representative B = p{B\) = p(-E>2)- The code 
B\ passes the parent test and therefore there is a canonical permutation <j>i : Bi — > B such 
that the last two coordinates of B\ are special with respect to (j>\. We can take for a canonical 
permutation of B2 the permutation 4>2 = 0i7? -1 , since 0i7r — 1 (-B2) = <fti(Bi) = p{B\) = p{B?), 
where 7? G S 2r +2, 5f(i) = 7r(i) for i € {1, 2, ... , 2r}, 7r(2r + 1) = 2r + 1, 7r(2r + 2) = 2r + 2. 
Then 

{0 2 (2r + l),^ 2 (2r + 2)} = {0x(2r + 1), 0i(2r + 2)} = {2r + l,2r + 2}, 

hence the last two coordinates of B2 are special with respect to the canonical permutation 
4>2- It turns out that the code B2 also passes the parent test. □ 

Theorem 9 If the set U s consists of all inequivalent binary self-dual [2s, s] codes, then the 
set U s +i obtained by the algorithm presented in Tabled consists of all inequivalent self-dual 
[2s + 2,s + l] codes, s > 1. 

Proof. We must show that the set U s +i filled out in Procedure Augmentation, consists 
only of inequivalent codes, and any binary self-dual [2s + 2, s + 1] code is equivalent to a code 
in the set U s+ \. 

Suppose that the codes B\, B2 € U s+ i are equivalent. Since these two codes have passed 
the parent test, their parents are also equivalent according to Lemma[71 These parents are self- 
dual codes from the set U s which consists only in inequivalent codes. We have a contradiction 
here and therefore the codes B\ and B2 cannot be equivalent. It follows that J7 s +i consists 
of inequivalent codes. 

Take now a binary self-dual code C of dimension s + 1 with a canonical representa- 
tive B. Hence B is equivalent to C and B passes the parent test. Since U s consists of all 
inequivalent self-dual codes of dimension s, the parent of B is equivalent to a code A G U s . 
According to Lemma [8] there is a child Ba of A, equivalent to B, such that Ba passes the 
parent test. Since the codes B and Ba are equivalent, so are the codes C and Ba- In this 
way we find a code in U s +i which is equivalent to C. □ 

Applying the algorithm recursively, we have the following corollary. 
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Corollary 10 If the set U r consists of all inequivalent binary self-dual [2r, r] codes, then the 
algorithm presented in Tabled generates all inequivalent self-dual [2k, k] codes, r < k. 

We can partition the set U r of all inequivalent binary self-dual [2r, r] codes into disjoint 
subsets U r i, U r 2, . . . , U rs and apply the algorithm to these subsets independently. Denote by 
Ufa the set of the inequivalent self-dual [2k, k] codes, obtained from C7 r j, i = 1,2, ... ,s, via 
the described algorithm. Following the algorithm and the theorems, we have 

Corollary 11 The union Uki U Uk2 U • • • U Uk s consists of all inequivalent binary self-dual 
[2k, k] codes, and U^i D E/jy = for i ^ j. 

This corollary shows that we can divide the computations into parts that need no 
mutual communication. 

The most difficult part in the algorithm is the calculation of the automorphism group 
and the canonical permutation of a given code. That's why we try to avoid this part by using 
invariants. Actually, we can do the parent test without knowing the automorphism group of 
the code, using only appropriate invariants, which is much faster. 

Definition 3 Let N = {1, 2, . . . , n} be the set of the coordinates of the code C. An invariant 
of the coordinates of C is a function f : N Z such that if i and j are in the same orbit 
with respect to Aut(C) then f(i) = /(j). 

The code C and the invariant / define a partition ir = {Ai, A2, ■ ■ ■ , iVj} of the coordinate 
set A, such that Aj C\Nj = for i ^ j, N = Ni U U • • • U N[, and two coordinates i,j are 
in the same subset of N <^=^> f(i) = f(j). So the subsets Ni are unions of orbits, therefore 
we call them pseudo-orbits. We can use the fact that if we take two coordinates from two 
different subsets, for example s £ Ni and t G Nj, Ni fl Nj = 0, they belong to different orbits 
under the action of Aut(C) on the coordinate set N. Moreover, using an invariant /, we can 
define a new canonical representative and a new canonical permutation of C. 

Firstly, we introduce an ordering of the pseudo-orbits of C. We say that Ni -< Nj for 
* + j, if: (1) | Ai| < \Nj\, or (2) |JV f | = [Nj\ and f(s) < f(t) for s G N h t G Nj. Then we 
redefine the canonical representative of one equivalence class in the following way: 

1. If the smallest pseudo-orbit consists of only one coordinate, we take for a representative 
a code in the equivalence class for which this coordinate is the last one, and the n — 1-th 
coordinate is from the second smallest pseudo-orbit. The last two coordinates of the 
code C have been added according the construction method described in the previous 
section. If none of these two coordinates belongs to the smallest pseudo-orbit, C does 
not pass the parent test. If one of them belongs to this orbit, we check whether the 
other one belongs to the second smallest pseudo-orbit. If no, C does not pass the parent 
test, but if yes, we should find the canonical representative to be sure whether C passes 
the test. 
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2. If the smallest pseudo-orbit N s contains three or more coordinates, we again take for a 
representative a code in the equivalence class for which the coordinates from the smallest 
pseudo-orbit are the last coordinates. According to Lemma [2] not all coordinates in N s 
are equal, so we can take two different coordinates in the end. 

3. The smallest pseudo-orbit N s consists of exactly two coordinates. If these two coordi- 
nates are different, we can take them for the last two coordinates. But if they are equal, 
we are looking for the next pseudo-orbit. There is only one self-dual code such that its 
coordinates can be partitioned in subsets of two elements such that both coordinates 
in each subset are equal, and this is the code n = 2k. But the automorphism group 
of this code is • Sk and it acts transitively on the coordinates. 

The complexity of the algorithm mainly depends on two of its steps. The first one is 
"find the set Child(A) of all inequivalent children of A\ For this step, we have as input a 
set L of generating elements of the automorphism group of the code A, and the set D of all 
binary odd- weight vectors with length dim A The algorithm splits the set of these vectors 
into orbits under the action of the group Aut(A). Any orbit defines a child of A, as different 
orbits give inequivalent children. This step is computationally cheap. The complexity is 
linear with respect to the product \L\ ■ \D\, \D\ = 2 dimA ~ 1 (for details see [23J). 

In the step " if B passes the parent test" , using a given generator matrix of the code B 
we have to calculate invariants, and in some cases also canonical form and the automorphism 
group Aut(-B). Finding a canonical form and the automorphism group is necessary when the 
used invariants are not enough to prove whether the code B pass or not the parent test. If the 
code B passes the parents test, the algorithm needs a set of generators of Aut(-B) for the next 
step (finding the children). For this step, we have to generate the set M w of all codewords of 
weight < w. The complexity here is 0(^™JJ^ ( dimB )). To determine the invariants, we use 
the set Md, where d is the minimum weight of the code B. We use mostly the invariants f\ 
and /2 defined as follows: 

/i:{l,2,...,n}->{0,l}, / 2 :{l,2,...,n}->Z, 

where n is the length of B. Moreover, if s = YlveM v i then fi(i) = 1 if and only if s is 
odd, and faii) = s, I < i < n, v = (t>i,t>2, • • • ,v n ) G To calculate the values of /i, the 
algorithm needs only \M^\ operations in the case of bitwise presentation of the codewords. 
For fi the algorithm uses 0{n\Md\) operations. We use also a vector valued invariant / 
such that f(i) = /bW)- If the code is not rejected with the invariants, the algorithm 

generates the smallest set M w of rank dimB. The corresponding to this set binary matrix is 
the input in the algorithm for a canonical form. For the complexity of this type of algorithms 
see (HE]. 

The described algorithm is implemented in the program Gen-self-dual-bin of the 
package Self-dual-bin written in Borland Delphi 6.0. 
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4 The Results 



To find all inequivalent self-dual codes of length 38, we begin from the set U\q of all inequiv- 
alent self-dual codes of length 32 and dimension 16. We partition U\q into three subsets and 
run them on three cores using a PC Intel i5 4 core processor. The number of all codes con- 
sidered in the program (these are the inequivalent children for each code in U\q U Un U U\%) 
is 2,338,260,952. For 151,016,675 of them, a canonical form is computed. The number of the 
obtained inequivalent codes in the set U\§ is 38,682,183. The calculations took about four 
days. 

Generator matrices of all inequivalent self-dual codes of length 38 are saved in three files 
in the form used in [6]. A compressed version of these files will be available on the web-page 
http://www.moi.math.bas.bg/~iliya/. A file with additional information which contains 
the number of codewords of weights 2, 4, 6, 8, and the order of the automorphism group for 
each code will be also available. This information allows to compute the expressions in the 
mass formulas. 

The possible weight enumerators of the self-dual codes of length 38 are given by the 
following formula 

W(y) =1 + ay 2 + (13a + (3)y 4 + (57 + 76a + 7(3 - 7 - 4% 6 + (228 + 260a + 17/3 - 7 + 28% 
+ (560a + 7(3 + 67 - 1365 + 1520)y 10 + (728a - 43/3 + 67 + 248<5 + 10032)y 12 
+ (364a - 77/3 - 157 + 1005 + 37620)y 14 + (85614 - 572a - 11(3 - 157 - 764% 16 
+ (127072 - 1430a + 99/3 + 2O7 + 528% 18 + • • • + y 38 , 

where a, (3, 7, S are integers. The numbers of inequivalent codes, the numbers of different 
weight enumerators and the numbers of different orders of the automorphism groups for each 
minimum weight d are given in Table [3j 



Table 3: Numbers of inequivalent codes of length 38 



d 


2 


4 


6 


8 


J) codes 


519492 


27463982 


10695965 


2744 


(t weight enumerators 


3504 


7176 


88 


2 


ft orders of Aut(C) 


799 


764 


75 


18 



The smallest order j)Aut s and the largest order jJAut; among the automorphism groups 
are listed in Table U] for each minimum weight d. 

We give one more table with some additional results about the automorphisms of the 
codes. We list in Table [6] the number of self-dual [38, 19, d] codes C such that p k divides the 
order of Aut(C) where p is a prime and A: is a positive integer. 
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Table 4: Orders of the automorphism groups 



d 



2 4 6 8 



JtAut. 



2 4 11 



tfAut, 



2 19 -19! 2 1A ■ 3 • 7 • 16! 1032192 504 



For the verification of our results, we use the mass formula (fTJ) and also the following 
corollary from one lemma of Thompson [24J : 

Theorem 12 [11] Let n and d be even positive integers and let U be a family of inequivalent 
self-dual codes of length n and minimum weight at most d. Then U is a complete set of 
representatives for equivalence classes of self-dual codes of length n and minimum weight at 
most d if and only if 



For the constructed codes we obtain the same values of the left and the right expressions 
in the formula ([3]) for n = 38, namely: 

(d=2) 19137697424578816915816164139573797711865999715625; 

(d=4) 2009458229580775776160697234655248759745929970140625; 

(d=6) 75153737786321014028410076576106303614497780883259375; 

(d=8) 1331294783643400819931835642205311664028246404217737500; 

(all) 27222898185745116523209337325140537285726884375 (formula ©). 

Remark 2 To calculate the sums in the mass-formulas ([I]) and ([3]), we use the package 
LONGNUM of S. Kapralov for calculations with large integers |14j . 

5 Conclusion 

The generation of all inequivalent binary self-dual codes of length n > 38, using only standard 
computer algebra systems, seems to be infeasible. That is why we use special algorithmic 
techniques to surmount difficulties and to classify codes even with PC's. In this work we 
describe the classification of the self-dual codes of length 38. 




(3) 
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Denote by SD k the number of all inequivalent self-dual codes of dimension k. Obviously, 



SDk J2k)\ • 

Consider the sequence a k = SD k (2k)\/ IltlC 2 * + X )' k = 2 > • • • • Looking at the already 
known classifications of binary self-dual codes, we can calculate the values of a k for k < 19. 
We list the integer part of these values in Table [5j Moreover, if we count the number of 
the inequivalent self-dual codes with a trivial automorphism group, we see that there are no 
such codes for length n < 32, but for the larger lengths this number increases very fast. For 
example, more than a quarter of all codes of length 38, namely 10140257 inequivalent codes, 
have a trivial automorphism group. So we have the following conjecture. 

Conjecture 1 The sequence {a k , k = 10, 11, 12, ... } is decreasing. 



Table 5: Values of [a k \ for k < 19 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


[ak\ 


2 


8 


48 


597 


3162 


18974 


70836 


230631 


353061 


464937 


k 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


L a fcJ 


327440 


194067 


57659 


13482 


2004 


273 


34 


7 


2 


? 



Some other problems we are going to attack are: 



• Classification of the binary self-dual codes of length 40. The doubly-even self-dual codes 
have been classified [2]. We have already started with the generation of the optimal 
[40, 20, 8] codes and the results will be ready soon. We have a lower bound on the 
TT 1 - (2* + 1) 

number of codes by — — - > 4, 585, 657, 509. According to Conjecture 1 and 

Table O we have 



Conjecture 2 The number of the inequivalent binary self-dual codes of length 40 are 

TT 19 (2 i 4- 1) 

at most 2 • lll=n , '- < 9, 171, 315, 020. 

40! 

• We are going to apply the described algorithm for quaternary self-dual codes. 
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Table 6: Number of codes C for which p k divides the order of their automorphism groups 



p k \d 


2 


4 


6 


8 


2 


519492 


27463982 


557127 


453 


2 2 


478434 


27463982 


89141 


111 


2 3 


467048 


17031875 


23498 


34 


2 4 


397699 


16153273 


8201 


1 


2 5 


375614 


9310617 


3788 


_ 


2 6 


303984 


8307428 


1733 


_ 


2 7 


274730 


4686394 


846 


_ 


2 8 


212982 


3854920 


414 


_ 


2 9 


183525 


2208800 


220 


_ 


2 io 


140158 


1629411 


134 


_ 


2 11 


112968 


970908 


79 


_ 


2 12 


85784 


641657 


48 


_ 


2 13 


65569 


387900 


26 


- 


2 14 


48579 


237143 


14 




2 15 


35399 


140837 


8 




2 16 


25470 


83500 


5 




2 17 


17908 


49360 






2 18 


12729 


29175 






2 19 


8838 


17009 






2 20 


6161 


9646 






2 21 


4162 


5612 






2 22 


2774 


3228 






2 23 


1812 


1754 






2 24 


1177 


996 






2 25 


790 


505 






2 26 


551 


235 






2 27 


345 


91 






2 28 


226 


33 






2 29 


156 


6 







p k \d 


2 


4 


6 


8 


2 3U 


87 


2 


_ 


_ 


2 31 


47 


1 


_ 


_ 


2 32 


16 


_ 


_ 


_ 


2 33 


7 


_ 


_ 


_ 


2 34 


3 


_ 


_ 


_ 


2 35 


1 


_ 


_ 


_ 


3 


132329 


2743510 


3916 


_ 

85 


3 2 


45728 


424433 


185 


7 


3 3 


17286 


93437 


26 


1 


3 4 


6463 


22451 


3 


_ 


3 5 


1763 


4315 


_ 


_ 


3 6 


256 


366 


- 


- 


3 7 


51 


52 






3 8 


8 


5 






5 


4209 


16112 


9 




5 2 


449 


624 






5 3 


52 


40 






7 


4270 


21981 


14 


5 


7 2 


467 


771 






7 3 


83 


83 






7 4 


16 


13 






7 5 


1 


2 






11 


47 


29 






13 


17 


7 






17 


4 








19 


1 




2 


1 


23 


6 


4 






31 


2 


2 
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Abstract 

An efficient algorithm for classification of binary self-dual codes is presented. As an 
. application, a complete classification of the self-dual codes of length 38 is given. 

>. 

<^> Index Terms: Self-dual codes, Classification, Isomorph-free generation 

CO 

o\ 

^ ! 1 Introduction 

The self-dual codes form one of the important classes of linear codes because of their rich 
algebraic structure and their close connections with other combinatorial configurations like 
block designs, lattices, graphs, etc. 

The classification of self-dual codes began in the seventies in the work of Vera Pless 
[20] . where she classified the binary self-dual codes of length n < 20. The method used in the 
beginning remained essentially the same throughout the succeeding classifications. This is a 
recursive classification which proceeds from smaller to larger length and codes are classified 
up to equivalence. The process begins with the formula for the number of all self-dual codes 
of length n called a mass formula. The number of the self-dual binary codes of even length 

nis N(n) = nrii _1 ( 2i + 

Throughout this paper all codes are assumed to be binary. Two binary codes are 
called equivalent if one can be obtained from the other by a permutation of coordinates. 
The permutation a £ S n is an automorphism of the code C, if C = cr(C) and the set of all 
automorphisms of C forms a group called the automorphism group of C, which is denoted 
by Aut(C) in this paper. If C has length n, then the number of codes equivalent to C 
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is n!/|Aut(C)|. To classify self-dual codes of length n, it is necessary to find inequivalent 
self-dual codes C\, . . . , C r so that the following mass formula holds: 

In the survey [T2J Huffman summarized the classification of all binary self-dual codes 
of length n < 36. As the complete classifications for lengths 34 and 36 was not done at that 
time, we present here a new version of Table 1 from |12| . but instead of the question marks 
we put the correct numbers. Moreover, we extend the table with two more lengths - 38 and 
40. The number of all inequivalent singly-even (Type I) and doubly-even (Type II) codes 
is denoted by |jj and jjjj, respectively. In the table d max j (d max ,ii) is the biggest possible 
minimum distance for which Type I (Type II) codes with the given length exist, and $ m axj 
{^maxjl) is their number. 

The classification of the self-dual codes of length 34 was given in [3j. Using the self- 
dual [34, 17, 6] codes, Melchor and Gaborit classified the optimal self-dual codes of length 
36, namely the [36,18,8] codes [19]. Recently, Harada and Munemasa in [TT] completed the 
classification of the self-dual codes of length 36 and created a database of self-dual codes |10| , 
The doubly-even self-dual codes of length 40 were also classified by Betsumiya, Harada and 
Munemasa [2j. Moreover, using these codes, they classified the optimal self-dual [38,19,8] 
codes. A classification of extremal self-dual codes of length 38 was very recently obtained 
also in [1] by somewhat different techniques. 

Actually, the construction itself is easy, the equivalence test is the difficult part of 
the classification. Because of that, for larger lengths the recursive constructions are used 
preferably heuristic, for building examples for codes with some properties. There are also 
many partial classifications, namely classifications of self-dual codes with special properties, 
for example self-dual codes invariant under a given permutation, or self-dual codes connected 
with combinatorial designs with given parameters. These classifications are not recursive but 
they use codes with smaller lengths, that is why the full classification is very important in 
these cases, too. 

There exist several methods to construct self-dual codes of length n + 2 from self-dual 
codes of length n. In pQ, the authors describe three such methods recalling them the recursive 
construction, the building-up construction and the Harada-Munemasa construction. The first 
one (recursive construction) gives all inequivalent [n + 2, n/2 + 1, d + 2] codes starting from 
the self-dual [n,n/2,d] codes [19j. The other two constructions give all self-dual codes of 
length n + 2 from the self-dual codes of length n. Our construction is similar to the Harada- 
Munemasa construction but we offer a better way to deal with equivalence classes. 

In this paper, we present an algorithm for generating binary self-dual codes which gives 
as output exactly one representative of every equivalence class. To develop this algorithm, 
we use an approach introduced by Brendan McKay known as isomorph-free exhaustive gen- 
eration |18j . The constructive part of the algorithm is not different from the other recursive 
constructions for self-dual codes, but to take only one representative of any equivalence class, 
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we use a completely different manner. This approach changes extremely the speed of gener- 
ation of the inequivalent codes of a given length. Its special feature is that practically there 
is not equivalence test for the objects. 

As a result, the classification of all binary self-dual codes of length 38 is presented. The 
number of these codes is given in the following theorem and also listed in Table [TJ 

Theorem 1 There are 38 682 183 inequivalent self-dual codes of length 38. 

In Section [2] we present the theoretical foundations of the used construction method. 
The aim is to obtain all inequivalent self-dual codes of dimension k (length 2k) on the base 
of the inequivalent self-dual codes of dimension k — 1 (length 2k — 2). 

In Section E] we describe the used algorithm. Codes which are equivalent belong to the 
same equivalence class. Every code can serve as a representative for its equivalence class. We 
use the concept for a canonical representative, selected on the base of some specific conditions. 
This canonical representative is intended to make easily a distinction between the equivalence 
classes. 

In Section S] we give the results with some more information about the obtained codes. 

2 The Construction Method 

In this section we describe the theoretical foundations of the used construction method. The 
aim is to obtain all inequivalent self-dual codes of length n = 2k and dimension k on the base 
of the inequivalent self-dual codes of dimension k — 1. Proposition [3] gives us the possibility to 
develop a classification algorithm. Actually, this proposition presents a well known property 
of the self-dual codes but usually it is stated in the case d > 4 (see for example [H]). Here we 
give a variant for all d but for that we need the following lemma. We say that the coordinates 
i\, . . . ,i s , 1 < i\ < ■ ■ ■ < i s < n, of the code C are equal if Xi x = Xi 2 = ■ ■ ■ = Xi 3 for each 
codeword x = (x±,X2, ■ ■ ■ , x n ) G C. 

Lemma 2 If C is a binary self-dual code then C does not have more than two equal coordi- 
nates. In other words, if x^ = Xi 2 = • • • = Xi s for each codeword x = {x\,X2, ■ ■ ■ ,x n ) £ C 
then s < 2. 

Proof. If n = 2 then C = i% = {00, 11} and C has only two coordinates. Let n > 4. Since 
the dimension of C is n/2 > 2, not all coordinates are equal and s < n. Suppose that s > 3, 
for all codewords x £ C. In such a case the vector y = (yi, . . . , y n ) <G F2 
of weight 2 with y^ = y± 2 = 1 will belong to C 1 - = C. In this way we obtain a codeword for 
which yi x = yi 2 = 1 but yi A = - a contradiction. Hence s < 2. □ 

Any self-dual [n, n/2, 2] code for n > 2 is decomposable as i% ® C n _2 where C„_2 is a 
self-dual code of length n — 2. Furthermore, the number of the self-dual [n, n/2, 2] codes is 
equal to the number of all self-dual [n — 2, n/2 — 1] codes (n > 2). But even if a self-dual code 
of length n > 2 has minimum distance d = 2, it has two coordinates which are not equal. 
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Table 1: Binary self-dual codes of length n < 40 |12| 



n 


h 


hi 


1 1 I < > rt * ± 


LItd n y T 


(Jjt~r~i n T T T 
1 1 bLLJL * J J. 




Reference 


2 


l 




2 


1 






[20] 


4 


l 




2 


1 






[2D] 


6 


l 




2 


1 






[20] 


8 


l 


1 


2 


1 


4 


1 


[20] 


10 


2 




2 


2 






EQ] 

\ i 


12 


3 




4 


1 






PQ] 


14 


4 




4 


1 






[20J 


16 


5 


2 


4 


1 


4 


2 


m 


18 


9 




4 


2 






[20] 


20 


16 




4 


7 






[20] 


22 


25 




6 


1 






[22] 


24 


46 


9 


6 


1 


8 


1 


[22] 


26 


103 




6 


1 






0E1EI] 


28 


261 




6 


3 






0[H1[2I] 


30 


731 




6 


13 






0[H1[2I] 


32 


3210 


85 


8 


3 


8 


5 


[ami 


34 


24147 




6 


938 






[2] 


36 


519492 




8 


41 






mm 


38 


38682183 




8 


2744 






this paper, [T| (2j 


40 


? 


94343 


8 


? 


8 


16470 


[2J 



Proposition 3 Xei C be a binary self-dual [n = 2k > 2, k, d] code and Co = {(^i, • • • , x n ) G 
C, x n -\ = Xn}. If the last two coordinates of C are not equal and C\ is the punctured code 
of Cq on the coordinate set {n — l,n} then C\ is a self-dual [n — 2, k — 1, d\ > d — 2] code. 

Proof. Let n > 2 and the last two coordinates of the code C are not equal. This means 
that there is a codeword y € C such that y n -i ^ Un- Therefore Co is a subcode of C with 
dimension k — 1, and (00... 011) C. The punctured code of Co on the coordinate set 
{n — 1, n} is 

Cl = {(XI,X 2 , ■ ■ ■ ,X n - 2 ) | X = (X1,X 2 , • • • ,X n -2,X n -l,X n ) G C }. 

Since (00... 011) C, dimCi = dim Co = k — 1. Moreover, if x = (x\,X2, ■ ■ ■ ,x n -2,a,a), 

V = (yi,V2, ■■■ , Vn-2, b, b) G C then + x 2 y2 H + x n _ 2 y„-2 + oh + ab = 0. Hence 

xiyi + x 2 y 2 H h x n -2Vn~2 = for any two codewords (xi, . . . , x n _ 2 ), (y l5 . . . , y n _ 2 ) G Ci. 

It follows that the code C\ is self-orthogonal, and since its dimension is a half of its length, 
it is self-dual. □ 
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For our construction, we use the code G\ of dimension k — 1 to obtain a self-dual code 
of dimension k and length 2k. We describe this in the following corollary. Let C\ be a binary 
self-dual [2k — 2, k — 1] code with a generator matrix Gi. Let us extend Gi with two equal 
columns (a\,a2, ■ ■ ■ , afc_i) T and consider the code Go generated by the matrix 



We can choose the vector (a%, a®, . . . , afc-i) such that (11... 11) G Go- Obviously, Go is a 
self-orthogonal [n = 2/c,/c — 1] code and its dual code is Cq = (Go, (00 . . . Oil), x) where 
x G Cq \ (G , (00 . . . Oil)) and wt(x) is even. Consider the last two coordinates of the vector 
x. If x n —i = x n then (x\, X2, ■ ■ ■ , x n —2 ) G Ci = Gi and therefore x G (G , (00 . . . Oil)). 
Hence x n ^\ ^ x n . The next corollary follows immediately. 

Corollary 4 The code C = {Cq,x) = Go U (x + Go) is a binary self-dual [2k, k] code. 

We will call the code G a child type code (CTC) of Gi, and C\ - a parent type code 
(PTC) of G. Actually, we have two possible self-dual codes obtained in this way from C\, 
namely (Go, x) and (Go, (00 . . . 011) + x), but both codes are equivalent. 

Remark 1 To describe the search tree for our algorithm, we use the terms CTC and PTC. 
In the literature on algorithms |15t [18] the terms are parent and child, but as these terms 
are used in a different way in some papers and chapters devoted to self-dual codes [13], we 
decided to set child type code and parent type code. 

In this way all self-dual codes of length n = 2k can be constructed on the base only on 
the inequivalent self-dual codes of length n — 2. Indeed, if we take two equivalent self-dual 
codes C\ = C[ and use the same vector (ai, a%, ... , afc-i), we can obtain equivalent self-dual 
codes of length 2k via the described construction. Let G\ be a generator matrix of C\ and 
P be a permutation matrix such that G\P generates the code C[. Then the matrix 



generates a self-dual code equivalent to G, namely the code C'. Furthermore, if it G SVi-2 is 
the permutation corresponding to the matrix P then tt G S n sends G to C' where Tr(i) = ir(i) 
for 1 < i < n — 2, and n(i) = i for i G {n — 1, n}. 

Let see now what happens if we take different vectors a and b of length k — 1 and use 
them in the described construction for the same self-dual [2k — 2, k — 1] code C\ with a gener- 
ator matrix G\. Consider the elements of the automorphism group Aut(Gi) as permutation 
matrices of order n — 2. To any permutation matrix P G Aut(Gi) we can correspond an 
invertible matrix Ap G GL(A: — 1,2) such that Gj = G\P = ApGi, since G^ is another gen- 
erator matrix of C\. In this way we obtain a homomorphism / : Aut(Gi) — > GL(/c — 1, 2). 
Consider the action of Im(f) on the set Fg -1 defined by A(x) = Ax T for every x G F^ 1 . 
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Theorem 5 [llj The matrices (G\ a T a T ) and (G\ b T b T ) generate equivalent codes if and 
only if the vectors a and b belong to the same orbit under the action of Im(f) on Fg -1 . 

Proof. Let the matrices [G\ a T a T ) and (G± b T b T ) generate the codes Co and C , respectively, 
and a T = A P b T , where P G Aut(Ci). Then 

T T\ ( P \ _ (n D „T „T\ _ ( a n a uT a uT\ _ a lr< u? uT\ 



(Gi a 1 a 1 )^ Q j 2 J = (GiP a 1 a 1 ) = {A P G X Apb 1 Apb 1 ) = Ap{G x b 1 b 

Since Ap{G\ b T b T ) is another generator matrix of the code C' , both codes are equivalent. 

Conversely, let Cq = C' Q . It turns out that there is a matrix B G GL(A;, 2) and an n x n 
permutation matrix P such that (G\ a T a T ) = B{G\ b T b T )P = BA P (Gi b T b T ). Hence 
BApGi = G\ and therefore BAp defines an automorphism of the code C\ = (G±). Since 
a T = BApb T , the vectors a and b belong to the same orbit under the action of Im(f) on 

F^r 1 . □ 

fx 1 \ 

Let now G = ( ^ ^ T ^ T J be a generator matrix of the self-dual [n, n/2, d] code C. 
If P € Aut(Ci) and y = xP then 



P \ _ / x 10\(P0\_fxP 10 
I 2 J ~\G 1 a T a T ){ I 2 ) ~ { G X P a T a T 



y 10 \_fl 0...0 \ f y 10 

A P d A P b T A P b T )- V T ^ P y V G x b T b T 

Hence the code C is equivalent to the code generated by the matrix G' = f ^ ^ T j , 

The following proposition reduces the number of the considered cases. It is a particular 
case of Theorem 1 from [161. 



Proposition 6 If C is a binary self-dual [n = 2k > 2, k, d] code then C is equivalent to a 
code with a generator matrix in the form 



G 



( x x ... 00 ... 1 \ 

X\ X\ 

4-i A \ \ 



and the matrix {Ik-\\A) generates a self-dual [n — 2,k — 1] code. 

Proof. If the last two coordinates of C are equal then according to Lemma [2] the first and 
the last coordinates are not equal and we can transpose them. So we can consider the case 
when the last two coordinates of C are not equal without loss of generality. 
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Take a generator matrix in systematic form for the code C\ obtained from C by the 
construction in Proposition [3l Then 



is a generator matrix of Cq in systematic form. Since Co is a self-orthogonal [2k, k — 1] code, 
Cq = (C , (00 . . . Oil), y) where y G C \ C . Hence y n _i ^ y n and suppose y n _i = 1, y n = 0. 

Consider the vector x = (xi,X2, ■ ■ ■ , Xfc_i, 00 . . . 010). Obviously, x 1 Co and hence 
x € Cj". But x C , x (00... 011) + C , x (00... 011) + y + C and therefore 
x G y + Co C C. It turns out that x G C and we can take a generator matrix for C in 
the needed form. □ 

If we start with a self-dual [2k — 2, k — 1] code C\ with a generator matrix in systematic 
form, for the construction in the last theorem we have to use a vector (xi,X2, . . . ,Xk-i) of 
odd weight. Then we are sure that the all-ones vector (11 ... 11) belongs to the subcode Co 
and the code C generated by the matrix ([2]) is a self-dual [2k, k] code. 

Theorem [5] and Proposition [6] are very important for our search. We use them to write 
a recursive algorithm which gives us all self-dual [n,n/2] codes starting from the code 12- 



In this section we present an algorithm for generating self-dual codes of a given length n. We 
discuss the background and give preliminary definitions and notations. 

Let f2fc be the set of all binary self-dual codes of dimension k (and length 2k). We 
consider the action of the symmetric group S^fc on the set f^, k = 1,2,... . This action 
induces an equivalence relation in as two codes Ci,C2 G £lk are equivalent (C\ = C2) if 
they belong to the same orbit. Hence the equivalence classes for the defined relation are the 
orbits with respect to the action of the symmetric group. According to Proposition El if we 
take two nonproportional coordinates in any code belonging to O^, k > 2, we can obtain a 
code from f2fc_i. Conversely, if we take all codes from £lk-i and extend them using Corollary 
[i] in all possible ways, we will obtain all codes from (7^. The construction of the self-dual 
codes of dimension k using the codes from fifc_i seems to be trivial, but actually this is a 
difficult problem and the question is how to find only the inequivalent codes. To do this, we 
develop a McKey type algorithm for isomorph-free generation |15} [TB] . 

The algorithm is an exhaustive search over the set of codes f2 = riiUr^U-'-Uf^ 
and this set is our search space. The generation process is described by a rooted tree (or 
forest). The nodes in the tree are objects from the search space fL From a self-dual code A 
of length 2k — 2 corresponding to the node A, we obtain self-dual codes of length 2k which 
are child type codes of A. To construct the child type codes, we use a generator matrix Ga 



\ 



A 




3 The Algorithm 
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of A in systematic form and the binary vectors from F^ -1 of odd weight. We denote the set 
of inequivalent child type codes by Child(A). The elements of Child(A) correspond to the 
nodes of the next level which are connected to A by edges. To find only the inequivalent 
child type codes, we use Theorem [5j Practically the rule A — > Child(A) which juxtapose 
child type codes to a code defines our search tree. The execution of the algorithm can be 
considered as walking the search tree and visiting all nodes through the edges. This can be 
done by a depth first search. We present a pseudocode of the algorithm in Table [2 

The following definitions and theorems help us to explain the algorithm and to prove 
its correctness. 

Definition 1 A canonical representative map for the action of the group S2k on the set Sif- 
ts a function p : Qk — > fl k that satisfies the following two properties: 

1. for all X £Q k it holds that p(X) = X, 

2. for all X,Y en k it holds that X = Y implies p(X) = p(Y). 

For a code CGflj, the code p(C) is the canonical form of C with respect to p. Analo- 
gously, C is in canonical form if p(C) = C. The code p{C) is the canonical representative of 
its equivalence class with respect to p. 

We can take for a canonical representative of one equivalence class a code which is 
more convenient for our purposes. Suppose that A is the canonical representative of this 
class with respect to a canonical representative map p. We can take B = A for a canonical 
representative for this class if we change the canonical representative map in the following 
way: p'{X) = B if X ^ B, p'{X) = p(X) if X ¥ B. According to Lemma [2] a self-dual code 
does not have more than two equal coordinates, hence if k > 2 we can take for a canonical 
representative of any equivalence class a code for which the last two coordinates are not equal. 

Example 1 If we order the codewords in any code lexicographically and then compare the 
codes according to a lexicographical ordering of the columns, we will have one biggest code 
in any equivalence class. We can take this code as a canonical representative of its class. 

Definition 2 Let C be a self-dual code of length n and p(C) = C be its canonical form. A 
permutation <pc E S n is called a canonical permutation of the code C, if 4>c(C) = p{C). 

For a fixed canonical representative map p, the canonical permutation of C depends on 
the automorphism group of the code and the permutations from the coset (/>cAut(C) can also 
be canonical permutations since 4>ccr(C) = 0c(o"(C)) = 4>c{C) = p(C) for any a E Aut(C). 

We obtain the automorphism group and a canonical form of a given code A using a 
modification of the algorithm presented in the paper |5]. This algorithm gives the order of 
the group, a set of generating elements, and a canonical permutation. 

The parent test can be defined in the following way. Suppose that A is a self-dual code of 
dimension k — 1 and B € Child(A). This means that B is obtained from A by the construction 
of Proposition El The parent test for B depends on the automorphism group of the code B. 
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Let ci and C2 be the coordinate positions of the code B for which 4>b(ci) = 2A;, 4>b{c2) = 2k — 1 
where 0b is the canonical permutation of B. The corresponding two columns of B are not 
equal, because the last two coordinates of the canonical representative are not equal. We call 
the coordinates c\ and C2 special for the code B with respect to the canonical permutation 
(pB- If there is an automorphism a of B such that {<t(ci), 17(02)} = {2& — 1,2/c} then the 
code B passes the parent test. In such case we can change the canonical permutation, taking 
4>b&~ 1 instead of (\>b- Then the last two coordinates of B are special with respect to the new 
canonical permutation. So a code passes the parent test, if there is a canonical permutation 
<p for this code such that the last two coordinates are special with respect to (p. 

Obviously, the canonical representative B of one equivalence class passes the parent 
test. If B\ = B also passes the parent test, then there is a permutation <j) : B — > B\ such that 
1 < (f>{i) < n — 2 for 1 < i < n — 2 and {<j)(n — 1), (j)(n)} = {n — 1, n}, where n is the length 
of the codes. This means that the parent type codes of these two codes are equivalent, too. 
So we have the following lemma. 

Lemma 7 If Bi and B2 are two equivalent self-dual [2k, k] codes which pass the parent test, 
their parent type codes are also equivalent. 

Lemma 8 Let A± and A2 be two equivalent self-dual codes of dimension r. Then for any 
child type code B\ of A\ which passes the parent test, there is a child type code B2 of A2, 
equivalent to B\, such that B2 also passes the parent test. 

Proof. Let G\ be a generator matrix of A\ in systematic form and B\ be the code obtained 
from A\ and the vector a = (a\, . . . , a r ) by the construction described in Corollary [3] and 
Proposition [6j Let B2 be the code generated by the matrix 7r(Gi) and the same vector a, 
where ir E #2r and tt(Ai) = A2- Obviously, 7r(Gi) generates the code tt(Ai) = A2. Moreover, 
according to Theorem [5j B2 is equivalent to all codes obtained by A2 and the vectors from 
the orbit with representative a under the action of Aut(^2) on Frj. Since the codes B\ and B2 
are equivalent, they have the same canonical representative B = p{B\) = p(B2). The code 
B\ passes the parent test and therefore there is a canonical permutation (f)\ : B\ — ^ B such 
that the last two coordinates of B\ are special with respect to <f>\. We can take for a canonical 
permutation of B2 the permutation §2 = ^it? -1 , since <\>yk~ x {B2) = ^i(-Bi) = p{Bi) = p{B2), 
where 7? G S 2r +2, = tt(>') for ie {1,2,..., 2r}, vr(2r + 1) = 2r + 1, 7? (2r + 2) = 2r + 2. 
Then 

{</> 2 (2r + l),<M2r + 2)} = {0i(2r + 1), 0i(2r + 2)} = {2r + l,2r + 2}, 

hence the last two coordinates of B2 are special with respect to the canonical permutation 
4>2- It turns out that the code B2 also passes the parent test. □ 

Theorem 9 // the set U s consists of all inequivalent binary self-dual [2s, s] codes, then the 
set U s +i obtained by the algorithm presented in Table\M consists of all inequivalent self-dual 
[2s + 2, s + 1] codes, s > 1. 
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Table 2: The main algorithm 



Procedure Augmentation^: binary self-dual code); 
begin 

If the dimension of A is equal to k then 
begin 

U k := U k U {A}; 

PRINT (A,tfAut(A)); 
end; 

If the dimension of A is less than k then 
begin 

find the set Child(A) of all inequivalent child type codes of A; 

( using already known Aut(-A);) 
For all codes B from the set Child(A) do the following: 
if B passes the parent test then Augmentation(i?); 

end; 

end; 

Procedure Main; 

Input: U r - nonempty set of binary self-dual [2r, r] codes; 

k - dimension of the output codes (k > r); 
Output: Uk - set of [2k, k] binary self-dual codes; 
begin 

Uk ■= (the empty set); 

for all codes A from U r do the following: 

begin 

find the automorphism group of A; 
Augmentation(A) ; 
end; 
end; 



Proof. We must show that the set U s +i filled out in Procedure Augmentation, consists 
only of inequivalent codes, and any binary self-dual [2s + 2, s + 1] code is equivalent to a code 
in the set U s +i. 

Suppose that the codes B\,Bi G U s +i are equivalent. Since these two codes have 
passed the parent test, their parent type codes are also equivalent according to Lemma \7\ 
These parent type codes are self-dual codes from the set U s which consists only in inequivalent 
codes. We have a contradiction here and therefore the codes B\ and B<i cannot be equivalent. 
It follows that U s+ \ consists of inequivalent codes. 
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Take now a binary self-dual code C of dimension s + 1 with a canonical representative 
B. Hence B is equivalent to C and i? passes the parent test. Since U s consists of all in- 
equivalent self-dual codes of dimension s, the parent type code of B is equivalent to a code 
A £ U s . According to Lemma [H there is a child type code Ba of A, equivalent to B, such 
that Ba passes the parent test. Since the codes B and Ba are equivalent, so are the codes C 
and Ba- In this way we find a code in U s +\ which is equivalent to C. □ 

Applying the algorithm recursively, we have the following corollary. 

Corollary 10 // the set U r consists of all inequivalent binary self-dual [2r, r] codes, then the 
algorithm presented in Tabled generates all inequivalent self-dual [2k, k] codes, r < k. 

We can partition the set U r of all inequivalent binary self-dual [2r, r] codes into disjoint 
subsets U r i, U r 2, • • • , U rs and apply the algorithm to these subsets independently. Denote by 
Uki the set of the inequivalent self-dual [2k, k] codes, obtained from U r i, i = 1,2, ... ,s, via 
the described algorithm. Following the algorithm and the theorems, we have 

Corollary 11 The union U^i U 11^2 U • • • U U^s consists of all inequivalent binary self-dual 
[2k, k] codes, and UkiCiUkj = for i ^ j. 

This corollary shows that we can divide the computations into parts that need no 
mutual communication. 

The most difficult part in the algorithm is the calculation of the automorphism group 
and the canonical permutation of a given code. That's why we try to avoid this part by using 
invariants. Actually, we can do the parent test without knowing the automorphism group of 
the code, using only appropriate invariants, which is much faster. 

Definition 3 Let N = {1, 2, . . . , n} be the set of the coordinates of the code C. An invariant 
of the coordinates of C is a function f : N Z such that if i and j are in the same orbit 
with respect to Aut(C) then f(i) = f(j). 

The code C and the invariant / define a partition it = {Aq, N2, ■ ■ ■ , Ni} of the coordinate 
set N, such that Ni PI Nj = for i 7^ j, N = N\ U N2 U • • • U Ni, and two coordinates i, j are 
in the same subset of f(i) = f(j). So the subsets N are unions of orbits, therefore 

we call them pseudo-orbits. We can use the fact that if we take two coordinates from two 
different subsets, for example s € N and t G Nj, N n Nj = 0, they belong to different orbits 
under the action of Aut(C) on the coordinate set N. Moreover, using an invariant /, we can 
define a new canonical representative and a new canonical permutation of C. 

Firstly, we introduce an ordering of the pseudo-orbits of C. We say that A, -< Nj for 
i / j, if: (1) |Aj| < \Nj\, or (2) |A,| = \Nj\ and f(s) < f(t) for s G N h t G Nj. Then we 
redefine the canonical representative of one equivalence class in the following way: 
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1. If the smallest pseudo-orbit consists of only one coordinate, we take for a representative 
a code in the equivalence class for which this coordinate is the last one, and the n — 1-th 
coordinate is from the second smallest pseudo-orbit. The last two coordinates of the 
code C have been added according the construction method described in the previous 
section. If none of these two coordinates belongs to the smallest pseudo-orbit, C does 
not pass the parent test. If one of them belongs to this orbit, we check whether the 
other one belongs to the second smallest pseudo-orbit. If no, C does not pass the parent 
test, but if yes, we should find the canonical representative to be sure whether C passes 
the test. 

2. If the smallest pseudo-orbit N s contains three or more coordinates, we again take for a 
representative a code in the equivalence class for which the coordinates from the smallest 
pseudo-orbit are the last coordinates. According to Lemma [2] not all coordinates in N s 
are equal, so we can take two different coordinates in the end. 

3. The smallest pseudo-orbit N s consists of exactly two coordinates. If these two coordi- 
nates are different, we can take them for the last two coordinates. But if they are equal, 
we are looking for the next pseudo-orbit. There is only one self-dual code such that its 
coordinates can be partitioned in subsets of two elements such that both coordinates 
in each subset are equal, and this is the code i|, n = 2k. But the automorphism group 
of this code is ■ Sk and it acts transitively on the coordinates. 

The complexity of the algorithm mainly depends on two of its steps. The first one is 
"find the set Child(A) of all inequivalent child type codes of A" . For this step, we have as 
input a set L of generating elements of the automorphism group of the code A, and the set 
D of all binary odd- weight vectors with length dim A. The algorithm splits the set of these 
vectors into orbits under the action of the group Aut(A). Any orbit defines a child type 
code of A, as different orbits give inequivalent child type codes. This step is computationally 
cheap. The complexity is linear with respect to the product \L\ ■ \D\, \D\ = 2 dimA ~ 1 (for 
details see [23]). 

In the step " if B passes the parent test" , using a given generator matrix of the code B 
we have to calculate invariants, and in some cases also canonical form and the automorphism 
group Aut(-B). Finding a canonical form and the automorphism group is necessary when the 
used invariants are not enough to prove whether the code B pass or not the parent test. If 
the code B passes the parents test, the algorithm needs a set of generators of Aut(S) for 
the next step (finding the child type codes). For this step, we have to generate the set M w 

of all codewords of weight < w. The complexity here is 0(Ya=1 C^'T^))- ^° determine the 
invariants, we use the set Mj, where d is the minimum weight of the code B. We use mostly 
the invariants A and A defined as follows: 

A :{l,2,...,n}^{0,l}, / 2 : {1,2,... ,n} -> Z, 

where n is the length of B. Moreover, if s = Y2veMd Vi then fi(i) = 1 if and only if s is 
odd, and A(*) = s,l<i<n, v = (vx,v%, • • • ,v n ) € F^- To calculate the values of fx, the 
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algorithm needs only \M^\ operations in the case of bitwise presentation of the codewords. 
For ji the algorithm uses 0{n\M ( i\) operations. We use also a vector valued invariant / 
such that f(i) = /bW)- If the code is not rejected with the invariants, the algorithm 

generates the smallest set M w of rank dim B. The corresponding to this set binary matrix is 
the input in the algorithm for a canonical form. For the complexity of this type of algorithms 
see E|. 

The described algorithm is implemented in the program Gen-self-dual-bin of the 
package Self-dual-bin written in Borland Delphi 6.0. 

4 The Results 

To find all inequivalent self-dual codes of length 38, we begin from the set U\q of all inequiva- 
lent self-dual codes of length 32 and dimension 16. We partition U\q into three subsets and run 
them on three cores using a PC Intel i5 4 core processor. The number of all codes considered 
in the program (these are the inequivalent child type codes for each code in U\§ U Un U U±g) 
is 2,338,260,952. For 151,016,675 of them, a canonical form is computed. The number of the 
obtained inequivalent codes in the set U\g is 38,682,183. The calculations took about four 
days. 

Generator matrices of all inequivalent self-dual codes of length 38 are saved in three 
files in the form used in [6]. A compressed version of these files is available on the web-page 
http://www.moi.math.bas.bg/~iliya/. A file with additional information which contains 
the number of codewords of weights 2, 4, 6, 8, and the order of the automorphism group for 
each code is also available. This information allows to compute the expressions in the mass 
formulas. 

The possible weight enumerators of the self-dual codes of length 38 are given by the 
following formula 

W(y) =1 + ay 2 + (13a + /3)y 4 + (57 + 76a + 7/3 - 7 - 4% 6 + (228 + 260a + 17/3 - 7 + 28% 
+ (560a + 7/3 + 67 - 1366 + 1520)y 10 + (728a - 43/3 + 67 + 248<5 + 10032)y 12 
+ (364a - 77/3 - 157 + 1005 + 37620)?/ 14 + (85614 - 572a - 11/3 - 157 - 764% 16 
+ (127072 - 1430a + 99/3 + 2O7 + 528% 18 + • • • + y 38 , 

where a, /3, 7, 6 are integers. The numbers of inequivalent codes, the numbers of different 
weight enumerators and the numbers of different orders of the automorphism groups for each 
minimum weight d are given in Table [3j 

The smallest order jJAut s and the largest order jJAutj among the automorphism groups 
are listed in Table [4] for each minimum weight d. 

We give one more table with some additional results about the automorphisms of the 
codes. We list in Table [6] the number of self-dual [38, 19, d] codes C such that p k divides the 
order of Aut(C) where p is a prime and A: is a positive integer. 
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Table 3: Numbers of inequivalent codes of length 38 



d 


2 


4 


6 


8 


jj codes 


519492 


27463982 


10695965 


2744 


(t weight enumerators 


3504 


7176 


88 


2 


ft orders of Aut(C) 


799 


764 


75 


18 



Table 4: Orders of the automorphism groups 



d 


2 4 


6 8 


tlAut s 


2 4 


1 1 


tfAut, 


2 19 • 19! 2 13 • 3 • 7 • 16! 


1032192 504 



For the verification of our results, we use the mass formula (pQ) and also the following 
corollary from one lemma of Thompson [23] : 

Theorem 12 [11] Let n and d be even positive integers and let U be a family of inequivalent 
self-dual codes of length n and minimum weight at most d. Then U is a complete set of 
representatives for equivalence classes of self-dual codes of length n and minimum weight at 
most d if and only if 

n/2-2 

n^+i). (3) 

i=i 

For the constructed codes we obtain the same values of the left and the right expressions 
in the formula ([3]) for n = 38, namely: 

(d=2) 19137697424578816915816164139573797711865999715625; 

(d=4) 2009458229580775776160697234655248759745929970140625; 

(d=6) 75153737786321014028410076576106303614497780883259375; 

(d=8) 1331294783643400819931835642205311664028246404217737500; 

(all) 27222898185745116523209337325140537285726884375 (formula 

Remark 2 To calculate the sums in the mass-formulas (pQ) and ([3]), we use the package 
LONGNUM of S. Kapralov for calculations with large integers |14j . 



^^M l{xGC|wt(x) = d}l = ( 



14 



5 Conclusion 



The generation of all inequivalent binary self-dual codes of length n > 38, using only standard 
computer algebra systems, seems to be infeasible. That is why we use special algorithmic 
techniques to surmount difficulties and to classify codes even with PC's. In this work we 
describe the classification of the self-dual codes of length 38. 

Denote by SD k the number of all inequivalent self-dual codes of dimension k. Obviously, 

SDk - — m — • 

Consider the sequence a k = SD k (2k)\/U^(2 i + 1), k = 1, 2, . . . . Looking at the already 
known classifications of binary self-dual codes, we can calculate the values of a k for k < 19. 
We list the integer part of these values in Table EJ Moreover, if we count the number of 
the inequivalent self-dual codes with a trivial automorphism group, we see that there are no 
such codes for length n < 32, but for the larger lengths this number increases very fast. For 
example, more than a quarter of all codes of length 38, namely 10140257 inequivalent codes, 
have a trivial automorphism group. So we have the following conjecture. 

Conjecture 1 The sequence {a k , k = 10, 11, 12, ... } is decreasing. 



Table 5: Values of \_a k \ for k < 19 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




2 


8 


48 


597 


3162 


18974 


70836 


230631 


353061 


464937 


k 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


[ak\ 


327440 


194067 


57659 


13482 


2004 


273 


34 


7 


2 


? 



The doubly-even self-dual codes of length 40 have been classified in [2] . A lower bound 

TT 1 - (2* + 1) 

on the number of all self-dual codes of that length is given by — % ~ l - > 4, 585, 657, 509. 

According to Conjecture 1 and Table [5j we have 

Conjecture 2 The number of the inequivalent binary self-dual codes of length 40 are at most 

TT 19 (2 i + 1) 
2 • iil=iy , ; < 9, 171, 315, 020. 
40! ' ' 
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Table 6: Number of codes C for which p k divides the order of their automorphism groups 



p k \d 


2 


4 


6 


8 


2 


519492 


27463982 


557127 


453 


2 2 


478434 


27463982 


89141 


111 


2 3 


467048 


17031875 


23498 


34 


2 4 


397699 


16153273 


8201 


1 


2 5 


375614 


9310617 


3788 


_ 


2 6 


303984 


8307428 


1733 


_ 


2 7 


274730 


4686394 


846 


_ 


2 8 


212982 


3854920 


414 


_ 


2 9 


183525 


2208800 


220 


_ 


2 io 


140158 


1629411 


134 


_ 


2 11 


112968 


970908 


79 


_ 


2 12 


85784 


641657 


48 


_ 


2 13 


65569 


387900 


26 


- 


2 14 


48579 


237143 


14 




2 15 


35399 


140837 


8 




2 16 


25470 


83500 


5 




2 17 


17908 


49360 






2 18 


12729 


29175 






2 19 


8838 


17009 






2 20 


6161 


9646 






2 21 


4162 


5612 






2 22 


2774 


3228 






2 23 


1812 


1754 






2 24 


1177 


996 






2 25 


790 


505 






2 26 


551 


235 






2 27 


345 


91 






2 28 


226 


33 






2 29 


156 


6 







p k \d 


2 


4 


6 


8 


2 3U 


87 


2 


_ 


_ 


2 31 


47 


1 


_ 


_ 


2 32 


16 


_ 


_ 


_ 


2 33 


7 


_ 


_ 


_ 


2 34 


3 


_ 


_ 


_ 


2 35 


1 


_ 


_ 


_ 


3 


132329 


2743510 


3916 


_ 

85 


3 2 


45728 


424433 


185 


7 


3 3 


17286 


93437 


26 


1 


3 4 


6463 


22451 


3 


_ 


3 5 


1763 


4315 


_ 


_ 


3 6 


256 


366 


- 


- 


3 7 


51 


52 






3 8 


8 


5 






5 


4209 


16112 


9 




5 2 


449 


624 






5 3 


52 


40 






7 


4270 


21981 


14 


5 


7 2 


467 


771 






7 3 


83 


83 






7 4 


16 


13 






7 5 


1 


2 






11 


47 


29 






13 


17 


7 






17 


4 








19 


1 




2 


1 


23 


6 


4 






31 


2 


2 
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